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INTRODUCTION 


This  is  the  third  in  a series  of  reports  [1-2]  prepared  in  support  of 
the  atmospheric  modeling  portion  of  the  Army  Materiel  Command’s  Target 
Signature  Analysis  (TSA)  program  (formerly  called  the  Ground  Target 
Signature  Program).  The  aim  of  the  atmospheric  modeling  project  is 
to  predict  absorption  and  scattering  effects  of  the  atmospheric  medium 
on  the  transmission  of  ground  target  electromagnetic  (EM)  signatures 
in  the  spectral  region  0.4-14  ym,  for  a wide  range  of  meteorological 
conditions.  Special  attention  is  being  given  to  low  visibility  atmos- 
pheric states  which  can  exist  under  battlefield  conditions. 

Two  important  related  aspects  of  the  program  are  the  atmospheric  ex- 
tinction nnd  radiative  transfer  models.  These  models  will  be  integrated 
into  system  performance  models  for  the  evaluation  of  terminal  homing, 
surveillance,  and  target  acquisition  systems.  Previous  reports  [1-2] 
have  described  atmospheric  transmission  at  1.06  ym  and  have  presented  a 
generalized  method  for  the  determination  of  EM  attenuation  due  to  the 
process  of  molecular  absorption.  This  report  will  describe  the  Atmos- 
pheric Sciences  Laboratory’s  (ASL)  single  scattering  model,  which  was 
developed  for  input  to  multiple  scattering  codes  for  the  determination 
of  EM  extinction  caused  by  the  atmospheric  aerosol.  These  multiple 
scattering  codes  will  be  described  in  a separate  report.  Sample 
applications  will  be  made  to  the  cases  of  haze  and  dust.  In  addition, 
the  role  that  the  ASL  single  scattering  model  plays  in  the  overall 
radiative  transfer  problem  will  be  described. 

In  the  remainder  of  this  paper,  in  addition  to  the  description  of  the 
ASL  single  scattering  model,  the  following  topics  are  considered: 
a discussion  of  the  physical  and  optical  characteristics  of  the  atmos- 
pheric aerosol  relevant  to  any  application  of  the  model;  the  versatility 
of  the  ASL  model  through  its  application  to  the  atmospheric  conditions 
of  haze  and  dust;  and  the  role  that  the  ASL  model  plays  in  the  overall 
radiative  transfer  problem.  The  theoretical  basis  of  the  ASL  scattering 
model  is  given  in  the  Appendix. 


THE  ASL  SCATTERING  MODEL 


In  the  spectral  region  of  interest  (0.4-14  ym)  the  scattering  of  electro- 
magnetic energy  (EM)  by  the  atmospheric  gases  is  well  known  and  will  not 
be  further  discussed.  The  scattering  coefficient  for  a Rayleigh  atmosphere 
can  be  easily  calculated  [3].  The  question  of  EM  extinction  by  aerosols 
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is  an  entirely  different  matter,  however,  because  of  the  aerosols’  complexity. 
The  ASL  single  scattering  model  is  a first  step  in  addressing  the  aerosol 
problem.  It  utilizes  the  Mie  scattering  theory  together  with  various 
size  distribution  options  to  compute  the  necessary  parameters  for  input 
to  multiple  scattering  codes  for  application  to  low  visibility  atmos- 
pheric conditions.  The  main  parameters  of  interest  computed  with  the 
ASL  model  are  the  extinction,  scattering,  and  absorption  cross  sections, 
the  single  scattering  albedo,  and  the  single  scattering  phase  function. 

In  the  context  of  this  report  the  atmospheric  aerosol  is  understood 
to  be  a colloidal  system  in  which  the  dispersed  phase  is  composed  of 
suspended  and  precipitating  solid  and  liquid  particles,  and  the  dispersive 
medium  is  air.  We  are  primarily  concerned  in  the  long  run  with  aerosols 
in  which  the  dispersed  phase  (solid  or  liquid  particles)  gives  rise  to 
low  visibilities  in  the  atmosphere;  i.e.,  when  the  visibility  is  11  km 
or  less . 

In  order  to  define  unambiguously  the  state  of  aerosols  in  the  atmos- 
phere, one  must  specify  the  following:  (1)  size  distribution,  (2) 

number  density,  (3)  composition,  (4)  shape,  and  (5)  spatial  and  temporal 
distribution. 

The  size  distribution  is  an  important  factor  and  defines  certain  types 
of  atmospheres.  The  size  distribution  for  hazes  is  quite  different  from 
that  for  fogs  and  clouds,  with  the  peak  of  the  distribution  lying  in  the 
region  of  smaller  sized  particles.  A variety  of  distribution  functions 
such  as  the  Junge,  modified  gamma,  and  the  log-normal  can  be  used  to 
model  a particular  state.  The  distribution  function  n depends  on  the 
size  of  the  particles  and  the  altitude;  i.e., 

n = n(z,  r) , (1) 

and  represents  the  particle  partial  number  density  per  unit  increment 
of  the  particle  radius  r at  an  altitude  z.  The  total  number  density  N 
at  altitude  z is  given  by 


N(z) 


r)  dr, 


r 

m 


(2) 


where  r and  rw  are  the  minimum  and  maximum  radii  of  the  distribution, 
m M 


Determining  the  composition  of  aerosols  is  a more  difficult  procedure. 
One  doesn’t  necessarily  need  the  exact  chemical  composition  but  rather 
the  optical  properties  of  the  aerosols,  such  as  the  real  and  imaginary 
parts  of  the  refractive  index.  Samples  of  the  atmosphere  can  be  taken 
from  various  geographical  areas  and  experiments  performed  to  determine 
the  bulk  complex  refractive  index.  It  is  especially  important  that  the 
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imaginary  part  of  the  refractive  index  be  determined,  since  results 
of  calculations  show  that  the  radiation  field  is  quite  sensitive  to 
that  component  [4]. 


The  shape  of  particulates  in  the  atmosphere  is  very  difficult  to 
estimate.  For  liquid  aerosols  one  can  usually  assume  the  shape  to 
be  spherical,  but  for  the  dry  dust  aerosols  the  shape  may  be  quite 
irregular.  How  the  shape  affects  the  optical  properties  of  a 
distribution  of  particles  is  not  well  known,  although  it  is  known 
that  the  scattering  properties  are  different  from  those  of  a spherical 
particle  [5].  Much  more  theoretical  and  experimental  research  on  this 
problem  is  needed.  In  this  report,  all  particles,  liquid  or  solid, 
are  assumed  to  be  spherical. 

s 

Finally,  the  spatial  and  temporal  distributions  of  aerosols  are  needed 
to  define  unambiguously  the  atmospheric  state.  It  is  usually  assumed 
that  the  vertical  distribution  of  particle  number  density  does  not 
vary  in  time  and  is  a monotonically  decreasing  function  of  altitude, 
and  sometimes  an  exponential  decrease  of  particle  concentration  in 
the  size  range  from  10“2  to  1 pm  up  to  a height  of  4-6  km  is  used;  i.e., 

N(z)  = N(o)e_z/H,  (3) 

where  N(z)  and  N(o)  are  the  particle  concentrations  at  heights  z and 
ground  level,  respectively.  H is  an  empirical  constant  that  varies  with 
location.  More  modeling  efforts  are  needed  to  define  the  vertical  dis- 
tribution and  its  change  in  time  for  various  hazes,  fogs,  clouds,  dust, 
smokes,  etc.  Most  radiative  transfer  models  assume  a horizontally 
homogeneous  atmosphere  - one  in  which  the  optical  properties  do  not  vary 
within  a horizontal  plane.  In  some  cases,  especially  in  a highly 
localized  military  setting,  this  is  probably  a poor  assumption. 


How  does  a knowledge  of  aerosols  aid  in  radiation  studies?  In  the 
following  ways:  In  order  to  determine  the  aerosol  optical  depth 

t between  points  z and  z , one  needs  the  volume  extinction  coefficient 


x 


z 

2 

Y (z)dz,  z <z<z 
ext  i 2 

z 

1 


(4) 


The  volume  extinction  coefficient  y (z)  in  turn  depends  on  the  total 

ext 

number  density  N(z)  and  the  total  extinction  cross  section  C : 


y (z)  = N (z)C 
Textv  ' v 7 ext 


(5) 
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where  the  effective  total  extinction  cross  section  C is  given  by 

ext 


C 


ext 


C 

ext 


(r) 


n(r)dr 


r 

m 
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(6) 


in  which  n(r)  is  the  size  distribution  function  for  spherical  particles 

of  identical  optical  properties,  C (r)  is  the  extinction  cross  section, 

e xt 

and  the  size  distribution  has  been  normalized  to  unity.  The  ASL 
scattering  model  will  allow  C ^ to  be  calculated  for  variable  complex 

refractive  index  m(A)  and  size  distribution  n(r).  Here  A is  the  wave- 
length of  the  incident  radiation.  The  model  can  also  be  used  to 
calculate  the  single  scattering  phase  functions  for  any  size  distribution 
and  complex  refractive  index. 

In  addition,  the  ASL  scattering  model  can  be  used  to  calculate  the 

scattering  cross  section  C and  therefore  determine  8 (z) > the 

sea  sea 

volume  scattering  coefficient.  Thus,  the  single-scattering  albedo 

a)  can  be  calculated: 
o 


0) 

o 


6 (z)  C C 

sea  sea  sea 


y (z)  C C + C , 

ext  ext  sea  abs 


(7) 


The  difference  of  the  extinction  and  scattering  cross  sections  is  the 

absorption  cross  section  C . which  can  be  used  to  calculate  the 

abs 

absorption  coefficient  a , 

abs 

In  order  to  do  more  than  a simple  parametric  study  of  the  natural 
radiation  field  in  the  atmosphere  by  arbitrarily  varying  the  optical 
properties,  one  needs  realistic  estimates  of  the  optical  depth  and 
the  scattering  phase  function.  These  in  turn  require  realistic 
descriptions  and  estimates  of  the  kinds  and  quantities  of  particulates 
in  the  atmosphere.  Since  the  natural  atmospheric  aerosol  consists  of 
both  liquid  drops  and  irregular  fragments  of  solid  materials,  the 
theoretical  prediction  of  the  scattering  properties  of  the  natural 
aerosol  has  a number  of  severe  limitations  which  stem  from  the  usual 
lack  of  information  concerning  the  complex  refractive  index  of  the 
particles,  their  internal  homogeneity  and  isotropy,  and  their  size, 
shape,  orientation  and  number  density.  In  view  of  these  uncertainties, 
it  is  normal  to  assume  that  all  aerosol  constituents  are  spherical  and 
isotropic,  and  that  size,  number  density,  and  bulk  optical  properties  may 
be  measured  or  assumed.  The  reason  for  this  is  that  rigorous  scattering 
theory  for  nonspherical  scatterers  is  extremely  complex.  Consequently, 
the  ASL  scattering  model  employs  the  well-known  Mie  theory  of  scattering 
by  spherical  particles  of  arbitrary  size.  Elements  of  Mie  theory 
necessary  for  an  understanding  of  the  numerical  algorithms  used  in 
the  ASL  model  are  discussed  in  the  Appendix. 
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Aerosols  are  generally  described  for  Mie  computations  by  their  complex 
refractive  index,  particle  size  distribution,  and  number  density.  A 
useful  approach  in  analyzing  the  impact  of  these  complex  and  highly 
variable  factors  on  the  transmission  of  EM  energy  is  to  adopt  typical 
or  average  values  for  these  parameters  of  the  dispersed  phase  for  a 
representative  set  of  meteorological  conditions  and  to  then  adjust 
these  values  as  required  for  other  conditions. 

We  have  adopted  five  different  options  for  describing  the  particle 
size  distribution:  an  arbitrary  one,  log-normal,  Junge  [6],  Deirmendjian 

model  C [7],  and  generalized  Khrgian  and  Mazin  [8,  9],  The  best  known 
distribution  function  - the  normal  or  Gaussian  - was  not  adopted  here 
because  it  is  symmetrical  and  is  calculated  under  the  assumption  that 
negative  particle  radii  can  occur. 

The  arbitrary  distribution  allows  the  input  of  tabulated  distributions. 
The  radius  r and  its  associated  frequency  of  occurrence  n(r)  are  simply 
read.  This  may  be  useful  in  comparing  results  between  experimentally 
measured  distributions  and  theoretical,  functional  ones. 

The  log-normal  distribution  which  has  been  adopted  does  not  admit 
negative  values  of  particle  radii,  but  it  does  show  the  skewness 
of  naturally  occurring  aerosol  particle  populations.  In  the  case 
of  the  log-normal  distribution,  it  is  x = lnr  rather  than  r which 
is  normally  distributed,  with  mean  value  x and  dispersion  a2  x. 

That  is, 


n (r)dr  = 


ax/2n 


exp 


(x-x): 

~2^ 


>dx, 


but  dx  = d(lnr)  = — , so  that 

r 


n (r)  = 


ro 


/2tt 


exp 


(x-x) 

2a2 

x 


(8) 


(9) 


where 


M 

x = j (lnr)n (r)dr 


m 


and 


M 

a2  = (lnr-x)/n(r)dr. 

X J 


m 


(10) 


(11) 


5 


The  mean  value  x is  usually  written  in  terms  of  the  geometric  mean 

value  r as  follows: 
g 


x 


1 

N 


N 


I lnri 


In 


i=i 


1/N 


lnr  . 
g 


(12) 


In  addition,  the  standard  deviation  ax  is  usually  written  in  terms  of 
the  geometric  mean  standard  deviation  Og  as  follows: 

a = lna  . (13) 

x g 

Hence,  one  can  calculate  the  frequency  of  occurrence  of  particles  with 
radius  r by  the  relation 


n(r) 


1 


exp 

rlna  /2tT 
g 


1 ln(r/r  ) ’ 

2 lna 

L g J 


2 


(14) 


This  distribution  has  been  found  to  describe  stratus  cloud  droplets, 
rock  dust  formed  by  mechanical  grinding,  mists  produced  by  a disk 
atomizer,  and  coalescent  aerosols  of  NH  Cl  and  H SO  formed  by  mixing 
gaseous  reagents  [10].  4 24 


The  Junge  distribution  is  given  by 


n(r)  = Cr  (15) 

where  C and  3 are  constants.  This  distribution  is  used  widely  to 
describe  natural  aerosols  such  as  mists  and  dusts,  and  is  generally 
accepted  as  valid  for  0.1  ym<r<15  ym.  For  smaller  particles,  a 
constant  frequency  of  occurrence  is  generally  used: 

n(r)  = C1.  (16) 

The  constant  3 in  Eq.  (15)  ranges  between  2 and  5 for  natural  aerosols, 
depending  upon  the  total  particle  concentration.  The  constant  C in 
Eq.  (15)  may  be  determined  by  integrating  Eq.  (15)  over  its  size  dis- 
tribution range  and  equating  the  result  to  experimentally  determined 
number  concentrations.  Work  done  at  White  Sands  Missile  Range,  NM,  by 
Henley  and  Hoidale  [11]  indicate  that  for  dust  C may  vary  from  0.3  to 
30  cm"^,  corresponding  to  low  and  high  dust  concentrations,  respectively. 
In  a similar  way  the  constant  C1  in  Eq.  (16)  can  also  be  determined. 
Henley  and  Hoidale  [11]  found  that  for  dust  C1  varies  from  10^  to  10 Jcm 
for  low  to  high  dust  concentrations.  In  the  specific  dust  distribution 
example  discussed  in  the  next  section,  we  have  chosen  to  use  the  dis- 
tributions Eq.  (15)  and  Eq.  (16)  with  3 = 4,  C = 30  cm-3,  and  C1  = 
105cra”3.  This  corresponds  to  a high  dust  concentration  and  therefore 
to  low  visibility. 
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The  Deirmendjian  model  C haze  is  a modified  Junge  distribution.  Since 
it  too  has  found  wide  application  to  continental  hazes,  we  have  included 
it  in  our  model.  The  actual  distribution  used  is 


n(r)  = 0 r<0.03pm 

n(r)  = 2.251  x 104  0 . 03pm<r<0 . lpm  (17) 

n (r ) = 2.251r“4  0.1pm<rl<*>. 

The  generalized  Khrgian  and  Mazin  distribution  [12]  was  first  proposed  to 
describe  cloud  droplets.  It  is  given  by 

n(r)  = araexp{-br^} » (18) 


where  the  quantities  a,  b,  a,  and  y are  positive  and  real  constants. 
These  four  constants  make  the  distribution  general  and  powerful  and 
allow  one  to  model  a wide  variety  of  diverse  aerosol  phenomena. 

For  a particular  choice  of  a and  y,  the  constants  a and  b can  be 
uniquely  determined  in  terms  of  the  total  number  of  particles  N per 
unit  volume  and  the  critical  or  mode  radius  rc.  For  example, 


a— i y y 

ar  (a-ybr  )exp(-br  ) . 


(19) 


Setting  this  equal  to  zero  at  r = r * where  the  concentration  is  at  a 
maximum,  one  obtains  c 


(20) 


Fur  the  more , 
get 


integrating  Eq.  (18)  over  the  entire  range  of  radii,  we 


N = a 


I a 

r exp 


-brY 


dr. 


o 

Thus  using  the  notation  T for  the  gamma  function  we  have 


(21) 


(orKD/y,  . 

a = Nyb  /t (^y=-) , <22) 

which  shows  that  for  a particular  choice  of  a and  y and  the  use  of  Eq . 
(20),  the  constant  "a"  is  found  from  the  knowledge  of  N.  In  particular, 
with  the  choice  of  a = 6,  y = 1,  rc  = 4 pm,  and  N = 100  cm"*3,  Deirmendjian 
[7]  obtains  the  cumulus  cloud  model 


n(r)  = 2.373r6e  1 * 5rcm  "pm  1 , 


(23) 
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and  with  a = 1,  y =1/2,  rc  = 0.05  ym,  and  N = 100  cm  3,  he  obtains  the 
coastal  haze  model 


n(r)  = 5.33  x 104re"8-  944v/r  . 


(24) 


SAMPLE  APPLICATIONS  TO  HAZE  AND  DUST 


In  order  to  illustrate  the  varied  conditions  to  which  the  ASL  scatter- 
ing model  may  be  applied,  three  representative  aerosols  were  chosen  for 
analysis:  (1)  an  inland  continental  haze,  (2)  a coastal  haze  and  (3) 

dust.  The  coastal  aerosol  is  represented  by  the  Deirmendjian  haze  M, 
which  in  turn  is  a particular  Khrgian  & Mazin  with  a = 1 and  y = 1/2. 

The  Mie  parameters  were  integrated  over  the  distribution  between 
rj  = 0.005  ym  and  r0  = 4.0  ym  in  Ar  = 0.01  ym  increments.  For 
illustration  purposes  only,  the  haze  particles  were  treated  as  pure  water 
droplets  and  the  refractive  index  of  pure  water  was  taken  from  Hale  and 
Querry  [13].  The  total  number  density  of  particles  used  in  these 
calculations  was  100  cm”3*  The  results  are  presented  in  Figures  1 and 
2. 


Figure  1 is  a plot  of  the  Mie  absorption,  scattering,  and  extinction 
cross  sections  versus  wavelength  for  the  spectral  region  0.4-14  ym. 

There  are  regions  of  absorption  centered  around  3 ym  and  6 ym,  and 
starting  at  10.5  ym  and  increasing  steadily  to  14  ym.  In  the  visible 
region,  the  greatest  loss  is  due  to  scattering.  However,  scattering 
is  also  predominant  between  the  narrow  absorption  regions  at  3 and 
6 ym.  Figure  2 shows  the  resulting  transmittance  over  a 1-km 
horizontal  path. 

The  inland  continental  aerosol  is  represented  by  the  model  C haze  of 
Deirmendjian,  which  is  a modified  Junge  distribution.  The  Mie  parameters 
were  integrated  between  r^  = 0.03  ym  and  r^  = 5 ym  in  Ar  = 0.01  ym  in- 
crements. The  refractive  index  of  dust  was  taken  from  Ivlev  and  Popova 
[14].  The  total  number  density  of  particles  used  in  these  calculations 
was  1. SyS'lO^cm”^  corresponding  to  a ground  visibility  of  5 km  [15]. 

The  results  are  presented  in  Figures  3 and  4. 

Figure  3 shows  the  Mie  absorption,  scattering,  and  extinction  cross 
sections  versus  wavelength  for  the  spectral  region  0.4-14  ym.  Regions 
of  absorption  exist  at  3,  7,  9,  and  11.5  ym.  Scattering  dominates 
below  6 ym,  especially  in  the  visible  region. 
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WAVELENGTH  (um) 


Figure  1.  Plot  of  the  Mie  absorption  (C  ) , scattering  (C  ),  and 

aos  sea 

extinction  (C  ) cross  sections  vs  wavelength  for  a coastal  haze. 


WAVELENGTH  (um ) 

Figure  2.  Plot  of  transmittance  in  a coastal  haze  over  a 1-km 
horizontal  path  at  sea  level  vs  wavelength,  particle  number  density 
100  cm 
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Figure  3.  Plot  of  the  Mie  absorption  (C  , ),  scattering  (C  ),  and 

abs  sea 

extinction  (C  t)  cross  sections  vs  wavelength  for  an  inland  haze. 


Figure  4.  Plot  of  transmittance  in  an  inland  haze  over  a 1-km 
horizontal  path  at  sea  level  vs  wavelength;  visibility  = 3 km. 
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In  Figure  4 we  have  plotted  the  resulting  transmittance  versus  wavelength 
for  a 1-km  path  length.  Through  the  visible  region,  the  transmittance 
increases  essentially  linearly  from  T = 0.33  at  0.4  ym  to  T = 0.56  at 
0.7  ym. 

Dust  was  chosen  as  representative  of  a dry,  localized  aerosol.  The 
distribution  used  for  the  giant  and  large  particles  is  a power  law 
with  6=4  and  C=30  cm-3,  and  for  the  smaller  particles  a constant  with 
C1=105cm“3.  However,  since  the  size  distributions  are  internally 
normalized  to  unity,  number  density  information  must  be  provided. 

The  number  concentration  of  dust  particles  was  related  to  visibility 
by  Henley  and  Hoidale  [11]  in  the  following  manner: 

N = • 105cm~3  , (25) 


where  V is  the  visibility  in  km  and  N the  total  number  density.  Values  used 
for  the  refractive  index  of  dust  were  obtained  from  Ivlev  and  Popova  [14] 
and  a visibility  of  1 km  was  used.  The  results  are  summarized  in 
Figures  5 and  6. 

Figure  5 is  a plot  of  the  Mie  cross  sections  versus  wavelength  for  the 
0.4-ym  to  12-ym  spectral  region.  Note  the  absorption  peaks  at  7,  9, 
and  11.5  ym  which  influence  the  total  extinction  cross  section.  As 
expected,  scattering  is  dominant  below  2 ym,  especially  in  the 
visible  region.  Figure  6 shows  the  transmittance  versus  wavelength 
over  the  spectral  region  for  a 1-km  path  length.  The  transmittance 
is  zero  through  the  visible  to  1.5  ym.  It  then  increases  steadily 
to  T = 0.23  at  6 ym.  Between  6 and  12  ym,  absorption  bands  reduce 
the  values  of  transmission  to  T = 0.04  at  7 ym,  T = 0.01  at  9 ym, 
and  T = 0.18  at  11.5  ym. 


ROLE  OF  ASL  MODEL  IN  RADIATIVE  TRANSFER  THEORY 


The  electromagnetic  radiation  field  in  the  visible  and  infrared  spectral 
regions  of  the  Earth’s  atmosphere  depends  upon  three  factors:  (1)  the 

natural  and/or  artificial  sources  of  radiation,  (2)  the  state  of  the 
atmosphere,  and  (3)  the  surface  conditions. 

During  the  day  the  most  important  source  of  natural  radiation  in  the 
visible  and  near-infrared  part  of  the  spectrum  is  the  sun.  In  fact,  more 
than  99%  of  the  solar  radiation  is  contained  within  the  spectral  band 
0. 2-4.0  ym.  Out  to  approximately  15  ym  the  main  source  of  natural 
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Figure  5.  Plot  of  the  Mie  absorption  (C  , ),  scattering  (C  ),  and 

abs  sea 

extinction  (C  cross  sections  vs  wavelength  for  a dust  haze. 


Figure  6.  Plot  of  transmittance  vs  wavelength  in  a dust  haze  for  a 
1-km  horizontal  path  at  sea  level;  visibility  = 1 km. 
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radiation  is  the  Earth Ts  surface.  About  30%  of  the  incoming  solar 
radiation  is  reflected  by  clouds,  dust,  and  the  Earth’s  surface,  whereas 
20%  is  absorbed  by  the  atmosphere.  The  remaining  50%  is  absorbed  by 
the  Earth.  The  consequent  emission  spectrum  of  the  Earth  is  that  of  a 
black  body  with  an  effective  temperature  of  255°K  with  a peak  in  the 
spectrum  at  12  ym. 

During  the  night  the  predominant  sources  of  natural  radiation  are  the 
moon,  stars,  airglow  in  the  visible  and  near-infrared  portions  of  the 
spectrum,  and  thermal  radiation  in  the  12 -ym  region.  The  intensity  of 
the  full  moon  is  about  10”6  that  of  the  sun  and  varies  in  a known  manner 
with  the  time  of  the  month.  Secondary  sources  of  radiation  are  those 
arising  from  the  scattering  or  absorption  of  radiation  by  the  atmos- 
phere and  the  ground.  This  includes  illumination  of  the  surface  by  light 
reflected  by  clouds  and  the  atmosphere.  Artificial  sources  of  radiation 
include  lasers,  searchlights,  flares,  and  thermal  radiation  from  fires, 
engines,  and  other  man-made  devices.  The  transfer  of  all  this  radiation 
through  the  atmosphere  is  described  by  radiative  transfer  models,  while 
the  relevant  parameters  used  by  these  models  are  described  by  scatterin 
models  such  as  the  ASL  model. 

In  the  last  15  years  many  different  radiation  models  have  been  developed 
for  calculating  the  radiation  field  existing  within  the  atmosphere. 

To  illustrate  the  role  that  the  ASL  scattering  model  plays  in  radiative 
transfer  theory,  described  below  is  a radiative  transfer  model  developed 
by  Turner  [16]  which  is  useful  for  determining  spectral  radiance  and 
irradiance  in  real  atmospheres. 

At  the  present  stage  of  development,  the  Turner  model  will  allow  one  to 
calculate  the  spectral  sky  radiance,  spectral  path  radiance,  total 
spectral  radiance,  the  downwelling  diffuse  spectral  irradiance,  the 
upwelling  diffuse  spectral  irradiance,  and  the  total  (direct  plus 
diffuse)  downwelling  spectral  irradiance  at  any  altitude,  view  angle, 
solar  angle,  and  for  any  value  of  albedo.  Also,  any  wavelength  X except 
those  in  the  gaseous  absorption  bands  can  be  used,  as  can  a variety  of 
atmospheric  states. 

The  model  is  designed  for  a plane-parallel  medium  with  a very  distant 
point  source  of  radiation  incident  upon  the  top  of  the  medium.  This 
point  source  is  usually  taken  to  be  the  sun,  but  it  could  just  as  well 
be  the  moon,  a star,  or  any  other  source  of  radiation  which  approximates 
a distant  point  source.  The  extraterrestrial  irradiance  values,  EQ(X), 
at  the  top  of  the  atmosphere  were  taken  from  the  NASA  standard  tables 
of  Thekaekara  [ 17 ] . 

All  calculations  are  spectral;  i.e.,  the  input  data  and  results  are 
given  at  specific  wavelengths,  not  over  bands.  In  the  visible  part 
of  the  spectrum  this  is  not  too  important  since  the  spectral  variations 
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of  the  basic  optical  parameters  are  small.  One  could  easily  simulate 
the  effects  of  a band  - for  instance,  from  0.55  to  0.60  ym  - by  taking 
an  average. 

Next,  one  must  provide  an  estimate  of  the  horizontal  visual  range  at  the 
surface.  However,  it  is  not  necessary  to  make  use  of  visual  range 
estimates  if  one  has  optical  depths.  The  optical  depth  can  be  computed 
with  the  use  of  the  ASL  model.  Optical  depth  is  the  more  fundamental 
quantity  and  is  preferred  over  visual  range  or  turbidity.  It  should 
be  noted  that  the  optical  depth  is  needed  for  all  the  wavelengths  and 
altitudes  which  are  desired. 

The  calculations  of  the  radiometric  quantities  radiance  and  irradiance 
are  definitely  a function  of  the  reflectance  or  albedo  of  the  under- 
lying surface.  Hence,  one  must  define  the  spectral  albedo  of  the  sur- 
face and  the  spectral  reflectance  of  the  target  under  observation.  The 
spectral  albedo  for  aerosol  particles  can  be  calculated  with  the  aid  of 
the  ASL  model.  At  the  present  time  only  Lambertian  surfaces  are  used  with 
the  Turner  model,  but  more  general  reflectances  are  being  considered. 

The  altitudes  of  interest  must  be  specified.  Any  altitude  from  0 to 
50  km  can  be  chosen,  realizing  that  the  optical  depth  for  all  practical 
purposes  is  zero  at  50  km,  the  assumed  top  of  the  atmosphere. 

Finally,  one  must  provide  values  of  the  solar  zenith  angle  0q,  the 
zenith  or  nadir  view  angle  0,  and  the  relative  azimuth  angle  <j>  between 
the  solar  plane  and  the  viewing  plane. 

The  major  advantages  of  this  model  over  many  others  are  the  following: 

(1)  since  the  solutions  for  the  basic  radiometric  parameters  are  closed- 
form  expressions,  the  computational  time  is  greatly  reduced;  and  (2)  the 
model  can  more  easily  be  adapted  for  a variety  of  complex  environmental 
situations,  such  as  non-Lambertian  targets  and  vertical  and  horizontal 
inhomogeneities . 

It  should  be  emphasized  that  the  use  of  visual  range  as  an  index  of 
atmospheric  state  is  just  a convenience.  If  one  has  available  the  more 
fundamental  values  of  optical  depth  versus  altitude  versus  wavelength 
for  a variety  of  atmospheric  conditions,  then  these  can  be  used  in  the 
model.  The  single-scattering  phase  functions  (see  Appendix) , which 
depend  upon  wavelength,  altitude,  and  atmospheric  condition,  are  also 
external  inputs  to  the  Turner  model.  Hence,  the  important  parameters 
of  albedo,  optical  depth,  and  phase  function  are  not  calculated  by  the 
radiative  transfer  model,  but  must  be  introduced  to  the  program  by 
using  either  a scattering  model  or  experimental  values. 
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CONCLUDING  REMARKS  AND  RECOMMENDATIONS 


ASL  has  the  capability  to  model  a variety  of  atmospheric  aerosol  states 
and  can  calculate  or  determine  the  three  critical  optical  parameters  - 
optical  depth,  single-scattering  albedo,  and  the  single-scattering 
phase  function.  In  order  to  insure  compatibility  with  the  input 
requirements  of  the  radiative  transfer  model,  the  following  model 
outline  is  used:  As  input,  the  type  of  aerosol  (haze,  fog,  dust, 

etc.)  is  specified  by  choosing  realistic  particle  size  distributions, 
complex  refractive  indices,  and  particle  number  densities.  Mie 
scattering  theory  is  used  to  calculate  the  polydisperse  cross  sections 
C as  a function  of  wavelength  and  particle  radius  for  complex  refractive 
index  m(A) . The  combination  of  total  number  density  with  total  cross 
sections  determines  the  attenuation  coefficients:  absorption  aa^s , 

scattering  £sca,  and  extinction  Yext*  These  coefficients  are  then 
used  to  determine  transmittance.  The  single-scattering  phase  functions 
are  calculated  directly  from  a knowledge  of  the  spectral  complex 
refractive  index  and  size  distribution.  Finally,  an  integration  of 
the  attenuation  coefficients  over  various  altitudes  gives  the  optical 
depth  as  a function  of  wavelength  A,  size  distribution  n,  visibility 
V,  and  altitude  z.  Also,  the  single-scattering  albedo  u)q  can  be 
calculated  from  the  coefficients.  Because  of  the  complexity  of  the 
theoretical  modeling,  it  is  recommended  that  a combined  theoretical- 
experimental  investigation  using  cloud  chambers  and  laser  scattering 
experiments  be  carried  out  to  determine  attenuation  coefficients  and 
phase  functions  for  various  atmospheric  states.  Work  should  continue 
to  relate  the  calculated  results  to  a convenient  set  of  available 
measurable  parameters  such  as  relative  humidity,  temperature,  visual 
range,  transmissivity,  turbidity,  contrast  transmittance,  etc.,  in 
order  to  define  the  atmospheric  state  in  an  unambiguous  way.  The 
main  task  in  this  phase  is  to  study  the  correlations  among  the 
measurable  parameters  P^,  P2,  P3,  ...  P-^.  This  parameter  set  can  then 
define  a unique  atmospheric  state  characterized  by  n2 , n^,  . . .n^) . 

As  an  example,  the  numbers  n^  could  be  the  following: 

n\  -*->  real  refractive  index 
n0  imaginary  refractive  index 
n3  size  distribution  parameter  #1 
++  size  distribution  parameter  #2 
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size  distribution  parameter  // 3 


n 


5 


particle  shape  factor  //I 


ni+l 


particle  shape  factor  #2 


nj  vertical  profile  index 


n 


j+1 


horizontal  inhomogeneity  index 


Many  nj_  values  are  required  to  define  a unique  atmospheric  condition. 
For  applications  to  realistic  environmental  situations,  however,  it  is 
advisable  to  combine  the  n^  values  into  a smaller  set  which  describes 
a class  of  states;  e.g.. 


n > n3>  n4  ni  ni+l * * * nk  •••nj  n j+1 * * * nN 


in  which  the  new  set  of  parameters  (A  > • • • A , . . . • • • A^)  can  more 
easily  be  identified  with  the  measurable  parameter  set  (P^ , P2,...P^). 
Thus,  for  a given  set  of  ’'measurements’1  or  readings  in  the  field, 
one  can  specify  a subset  of  possible  atmospheric  states,  and  in 
connection  with  a suitable  radiative  transfer  model,  analysis  will 
allow  the  radiation  field  to  be  determined. 

Naturally,  the  ASL  investigation  of  the  radiation  field  of  the  atmos- 
phere is  not  independent  of  other  efforts  by  other  laboratories.  We 
need  reflectance  models,  radiation  measurements,  etc.  Thus,  it  is 
recommended  that  all  models  being  developed  under  the  TSA  program  be 
interfaced  in  this  phase  of  the  work,  and  the  resulting  combined 
models  validated.  These  validation  studies  should  include  both  field 
and  laboratory  measurements. 
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APPENDIX  A 


MIE  THEORY 


Mie  theory  [18]  predicts  the  scattering  by  and  the  absorption  in  an 
isolated,  discrete,  homogeneous,  isotropic  sphere  of  diameter  D with  a 
known  complex  refractive  index  m = n-ik  relative  to  the  surrounding  medium 
and  illuminated  by  monochromatic  radiant  energy  with  wavelength  X in  the 
surrounding  medium.  The  theory  is  given  in  detail  in  standard  texts  [3,  4, 
19-21]  and  need  not  be  repeated  here.  Instead,  those  elements  of  theory 
needed  for  an  understanding  of  the  numerical  algorithms  used  in  the  ASL 
model  are  included. 
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Scatterers  attenuate  beams  of  radiant  energy  by  scattering  some  of  the 
energy  into  directions  other  than  the  incident  or  forward  direction  and 
by  absorbing  some  of  the  incident  energy  within  the  body  of  the  particle. 
The  combined  effect  of  pure  scattering  by  the  particle  and  true  absorption 
within  the  particle  is  termed  extinction.  The  amount  of  extinction, 
scattering,  and  absorption  by  a single  particle  is  given  in  terms  of 
corresponding  equivalent  blocking  areas  or  cross  sections,  Cexf. , Csca, 
and  Caks,  respectively.  These  cross  sections  depend  only  on  the  refractive 
index  of  the  particle  m = n-ik  and  the  size  parameter  a = 27Tr/A,  where 
r is  the  particle  radius. 

The  transmission,  T,  of  a cloud  of  particles  of  geometric  depth,  d, 
and  number  density,  N,  is  given  by 

T = e_T,  (A-l) 

with  the  optical  depth,  t,  given  by 

t = NdC  , (A- 2) 

ext 

and  the  cross  sections  related  by 

C = C + C , . (A- 3) 

ext  sea  abs 


When  the  particle  refractive  index  is  real  (k  = o) , the  absorption  in  the 
particle  vanishes;  i.e.,  Ca^g  = 0.  The  balance  between  loss  by  scattering 
and  loss  by  absorption  is  frequently  characterized  by  the  albedo  of 
single  scattering  0)o  > given  by 


C C 

sea  sea 


Uo  C ’ + C , C 

sea  abs  ext 


(A-4) 


A scatterer  with  wn  = 1 has  no  absorption  and  is  termed  a conservative 
scatterer.  The  albedo  gives  the  probability  that  a photon  encountering 
the  scatterer  will  be  scattered  into  some  direction  including  the  incident 
direction. 

Although  the  extinction  by  a cloud  of  particles  is  correctly  given 
by  Eqs . (A-l)  and  (A-2) , two  implicit  assumptions  may  lead  to  improper 
use  of  the  equations.  The  optical  depth  t in  Eq.  (A-2)  does  not  include 
losses  caused  by  absorption  in  the  medium  surrounding  the  particles. 
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This  assumption  obviously  breaks  down  at  wavelengths  for  which  atmos- 
pheric gases  absorb  appreciably.  The  second  assumption  is  that  scattered 
photons  never  return  to  the  incident  direction.  This  effect  becomes 
increasingly  important  as  optical  depths  exceed  t = 0.1. 

A final  caution  should  be  noted  in  regard  to  absorption  within  the 
particle.  Although  absorption  within  the  particle  is  correctly  determined 
by  the  wavelength-dependent  imaginary  part  k of  the  refractive  index  m, 

the  explicit  mechanism  which  causes  the  absorption  is  usually  not  spec-  4 

ified.  Texts  may  be  consulted  for  specific  examples  [22].  Usually  the 

absorption  is  joule  heating,  and  it  is  sometimes  necessary  to  account 

for  the  isotropic  black  body  radiation  emitted  by  the  scatterer  when 

its  temperature  rises  above  that  of  its  surroundings.  There  may  also 

be  circumstances  when  quantum  transitions  occur  in  the  scatterer  followed 

by  emission  at  or  near  the  same  wavelengths.  It  is  incumbent  on  the  user 

of  the  numerical  algorithms  presented  here  to  properly  include  these 

effects  since  they  are  not  automatically  accounted  for  in  these  algorithms. 

All  scattering  properties  of  spheres  are  computed  from  n and  k,  and 
through  the  use  of  the  induced  electric  and  magnetic  multipole  moments 
of  the  sphere  a^  and  b^,  respectively.  The  moments  are  given  by 


^(ma)  ^£(a)-m  ^(ma)  Y^Ca) 


(A-5) 


a£  y^Cma)  C^(a)-m  ^(ma)  ^(a) 


and 


m ^(ma)  Y£(a)  - ^(ma)  ^(a) 


(A-6) 


b£  m ^(ma)  ^(a)  “ Y£(ma)  5[(a) 


The  prime  denotes  differentiation  with  respect  to  the  argument. 
The  4^(z)  and  £^(z)  functions  are  Ricatti-Bessel  functions  of  the 
firsthand  third  kind,  respectively,  and  are  related  to  spherical 
Bessel  functions  j^Cz)  and  n^(z)  by 


y£(z)  = z j£(z) 


(A-7) 


and 


^(z)  = zj£(z)-i  zn£(z)  , 


(A-8) 


where 


(A-9) 


and 


1/2 

= (37)  N£+1/2^z->* 


(A-10) 
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The  function  J£+i/2(z)  is  t^le  integral  order  Bessel  function; 

the  function  ^£4-2/2^^  t^ie  ^alf  integral  order  Neuman  function. 

The  extinction  cross  section  is  computed  from 

00 

c . = Ai  E (24+1)  Re  (a.+  b.)  (A-ll) 

exc  2tt  £=1  * 

and  the  scattering  cross  section  from 

00 

^sca  ■ sf  l-l  (2I+1)  d»t|2+l1>tl2]-  (A-12) 

The  various  cross  sections  are  the  basic  quantities  used  in  scattering 
problems,  but  they  are  not  the  quantities  usually  computed  directly 
from  Mie  algorithms.  Instead,  it  is  more  convenient  to  compute 
dimensionless  efficiency  factors  Qe  and  Q , , which  depend  on  n,  k, 

and  a,  and  which  are  multiplied  by  the  geometrical  sphere  cross  section 
to  obtain  the  true  cross  section  =,rcrzQ_£. 

Although  the  cross  sections  account  for  the  energy  removed  from  the 
forward  beam,  they  do  not  give  any  information  about  where  the  scattered 
photons  go.  This  information  is  contained  in  scattering  amplitudes  and 
intensity  factors  which  relate  the  flux  density  scattered  through  an 
angle  0 relative  to  the  incident  flux  density.  There  are  two  amplitudes, 
S^G)  and  82(6),  and  intensity  factors,  i ^ ( 0 ) and  i2(0),  which  correspond 
to  light  respectively  polarized  perpendicular  and  parallel  to  the  plane 
of  scattering  defined  by  the  direction  of  incidence  and  the  direction  of 
scattering.  The  ratio  of  incident  to  scattered  density  at  radius  r from 
the  sphere  is  given  for  perpendicular  polarization  by 


ir  ^2ii  (0) 

Iro 


(A-13) 


and  for  parallel  polarization  by 

h m A2i2  (6) 

I£o  4ir2r? 


(A-14) 
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The  intensity  factors  are  related  to  the  scattering  amplitudes  by 

ij(0)  = | S2 (0)  | 2 

and 

i2(e)  = |s2(e) |2. 

The  amplitudes  come  from  the  multipole  moments  through 


(0) 


= L 


1=1 


s [a  it  (0)  + b x (6)] 

«,(«.+!)  i a a sl  1 


and 


S2(0)  ” \=1  '(Si)  tbJl1I£(0)  + a£T£(0)]. 


and  angular  factors  ^(0)  and  t^(0)  defined  in  terms  of  associated 
Legendre  functions: 


7T  (0)  = P 1 (cos  0)  /sin  0; 


Ve) 


dP*(cos  0) 
d0 


The  amplitudes  have  relative  phase  6 = argSj  - argS2. 

Alternative  expressions  frequently  used  are 

dP^(cos  0) 
d(cos  0) 

and 

dn£(0) 

r o(0)  = cos  0 ' 7rp(0)  - sin‘  0 * — — * 

* x d (cos  0) 

where 

1 a Si 

P £ (cos  0)  = (cos20  -1) 

dcos*-  0 


(A-15) 

(A-16) 

(A-17) 

(A-18) 

(A-19) 

(A-20) 

(A-21) 

(A-22) 

(A-2J) 
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These  functions  satisfy  the  following  recurrence  relations: 


Ve) 


cos  0 


<vi} 

(t-D 


Vi(e)  - 


£ 

£-1 


\-2(9) 


(A-24) 


and 


T £ (Q ) = cos  0 £tt^  ( e ) - 7T^2(e)]  “ (2£-l)sin207Tj^_1(0)  + T£_2(0),  (A-25) 


with 


7T0(e)  = 0, 

TO(0) 

tti  (e)  = 1, 

Tj  (0) 

IT  (0)  = 3 cos  0, 

t2(0) 

0, 


cos  0 , 

3 cos  2 0. 


The  scattering  cross  section  measures  the  ability  of  a particle  to 
scatter  light,  and  it  is  to  be  expected  that  Csca  is  obtained  from 
an  integral  over  the  scattering  intensity  factors.  Equation  (A-12) 
follows  from 

1 

C = -7—  i (i  (0)  + i (6))  dcose.  (A-26) 

sea  4tt  j 1 2 

-1 

Although  the  intensity  factors  themselves  may  be  used  in  scattering 
calculations,  they  are  primarily  suited  for  computing  flux  densities, 
and  it  is  frequently  more  convenient  to  measure  and  compute  scattered 
light  in  terms  of  radiances.  Radiances  do  not  have  the  1/r  depend- 
ence, and  it  is  therefore  unnecessary  to  know  the  distance  from  the 
scatterer  to  the  detector  if  the  detector  field  of  view  is  small 
and  is  filled  by  the  scattering  cloud.  The  phase  function  p(0) 
gives  a radiance  I scattered  into  the  0 direction  in  terms  of  the 
radiance  Iq  incident  on  the  particle. 

The  phase  function  is  dimensionless  and  is  defined  here  as 


p(0) 


A2 

2*Cext 


ij(e) 


(A-27) 
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The  normalized  phase  function  p ( 0 ) d^2/ 4tt  gives  the  probability  of  a 
photon  being  scattered  through  an  angle  0 into  an  element  of  solid 
angle  d ft  = d(f)dcos0.  The  integral  of  the  normalized  phase  function 
is  the  single  scattering  albedo  go  , which  gives  the  probability  that 
the  photon  is  scattered: 


GO 

o 


_1 

4tt 


27T 


/ I 


p (0)  d<J>dcos0 




4ttC 

ext 


i 


i (0)+i  (0)  dcos0 

1 2 J 


(A-28) 


or 


GO 

O 


C /C  . 
sea  ext 


(A-29) 


The  phase  function  contains  a sum  over  the  polarization  states  implicit 
in  the  ii  and  i2  intensity  factors,  and  is  thus  unsuitable  for  describing 
the  polarization  of  the  scattered  light.  Stokes  vectors  and  Mueller 
matrices  [3,  23]  are  a particularly  convenient  way  to  represent 
polarized  light  beams  and  their  interaction  with  scatterers  and  optical 
devices.  In  this  context  the  phase  function  is  replaced  by  a particular 
Mueller  matrix  - namely,  the  phase  matrix. 

Stokes  vectors  and  Mueller  matrices  can  be  variously  defined.  Here  we 
use  the  notation  of  Van  de  Hulst  [3].  Van  de  Hulst  writes  the  Stokes 
elements  as  time  averages  of  horizontal  and  vertical  electric  field 
components : 


_ - ie , -ik  z+iu)t 

e*  ■ aie  ie  ° ; 

(A-30) 

E - a e‘iE2e-lkoZ+1“t; 
r r 

(A-31) 

I - I.  + I = a2n  + a2  ; 
SL  r £ r 

(A-32) 

Q = I - I = a2n  + a2  ; 
x i r £ r 

(A-33) 

U = 2a  a cos6; 
l r 

(A-34) 

V = 2ana  sin6 ; 
l r 

(A-35) 

6 = e ^ - e . 

(A-36) 

24 


Here  r and  Z refer  to  the  electric  fields  perpendicular  and  parallel 
with  the  plane  of  scattering,  k is  the  propagation  constant,  u>  is 
the  angular  frequency,  (an , ar)  are  the  positive  amplitudes  of  (E^,  ® )» 
respectively,  and  (e^,  £ 2 ) are  the  corresponding  phases  of  the  waves. 
Furthermore,  the  electric  vector  is  given  by 


E = Re[E££  + Err],  (A-37) 

and  the  sense  is  chosen  so  that  t x Z is  in  the  direction  of  propagation. 

As  viewed  from  a detector,  these  elements  in  general  represent  an 
elliptically  polarized  beam  with  its  major  axis  at  an  angle  x with  the 
horizontal  reference  axis  and  with  ratio  of  minor  to  major  axis  tangent 
3.  The  corresponding  elements  of  the  Stokes  vector  are 


I = 

a2; 

(A-38) 

Q = 

a2cos23cos2x; 

(A-39) 

u = 

a^cos23sin2x ; 

(A-40) 

v = 

a2sin2B. 

(A-41) 

For  plane-polarized  light,  3=0  and  V vanishes.  3 is  +tt/4  for  right 
circular  and  -tt/4  for  left  circular  polarized  light.  The  tangent  of 
the  phase  difference  6 (Eq.  (A-36))  is 

tan  6 = tan  2B/sin2x*  (A-42) 

Stokes  vectors  add  linearly,  so  'that  a partially  polarized  beam  is 
represented  by  the  sum  of  an  unpolarized  beam  {l  0,  0,  0}  with  a 
completely  polarized  beam  {l  Q,  U,  V}.  This  gives  an  inequality 

I2  = (Iu  + Ip)2  > Q2  + U2  + V2,  (A-43) 


which  becomes  an  equality  when  Iu  = 0 and  the  beam  is  completely 
polarized.  The  plane  of  polarization  is  given  by  U/Q  = tan  2y,  and 
the  degree  of  polarization  is  [Q2  4*  U2  4* 
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Mueller  matrices  are  four-by-four  matrices  which  transform  an  incident 
Stokes  vector  into  an  outgoing  Stokes  vector: 

{I,  Q,  u,  V}  = M{I0,  Q0,  u0,  v0}.  (A-44) 

When  the  Stokes  vector  represents  a radiance  and  M represents  a Mie 
scatterer,  M is  the  phase  matrix 


(l/2)(i2+il) 


M = 


ttC 


ext 


(1/2)  (i^) 


(l/2)(i2-ii)  0 

(1/2) (i  +i  ) 0 


0 


0 


i cos6 


3 


0 


0 

-i^sin  6 


(A-/i  5) 


0 


0 i^sintf 


i cos  6 

3 


1/2 

in  which  13=  (i^  * i 2 ) , and  6 is  the  difference  in  phase  between 

amplitudes  S ^ (0 ) and  82(6) • The  corresponding  phase  matrix  for  the 
Stokes  vector  {I^5  lr,  U,  V}  is 


M1 


ttC 


ext 


0 

0 


0 

0 


0 0 


i cos6  -i  sin6 
3 3 


0 o 


i sin6  i cos6 

3 3 


(A-46) 


The  two  phase  matrices  and  all  corresponding  Mueller  matrices  are  related 
by  similarity  transforms 

M'  = T_1  Mr,  (A-47) 

M = TM'T-1,  (A-48) 

with  transform  matrices 


T 


1 1 

1 -1 

0 0 

0 0 


0 0 

I 

I 

0 0 


1 0 ' 

0 1 S 


(A-49) 
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and 


T 1 


1/2  1/2  0 0 
1/2  -1/2  0 0 

0 0 10 

0 0 0 1 


These  matrices  also  transform  the  Stokes  vectors  themselves: 


{I£,  Ir,  u,  V}  = T Hi,  q,  V,  U} ; 
{I,  Q,  U,  V}  = T{ 1^ , Iv,  U,  V}. 


(A-50) 


(A-51) 

(A-52) 


The  orientation  of  the  horizontal  or  reference  axis  implicit  in  a 
Stokes  vector  is  quite  arbitrary,  and  it  is  often  necessary  to  transform 
to  a new  reference  axis.  This  may  be  done  by  rotating  either  the  Stokes 
vector  or  by  applying  the  rotation  matrix  to  the  Mueller  matrix  in  a 
similarity  transformation.  The  rotation  matrix 


L(4>) 


10  0 0 

0 cos2(f>  sin2(J)  0 


0 -sin2(J)  cos2(J)  0 


0 0 


0 1 


(A-53) 


rotates  the  axes  clockwise  through  the  angle  cf> . The  inverse  transform 
L”*((j))  is  L(-(j>)  , so  that  Mueller  matrices  transform  as 


\ = L(-*)M0  L(*). 


(A-54) 


Although  it  is  unnecessary  for  single  scattering  calculations,  multiple 
scattering  calculations  are  usually  done  in  a system  of  coordinates 
embedded  in  an  atmosphere  rather  than  attached  to  each  individual 
scatterer  or  to  a particular  source-scatterer-detector  combination. 
Accordingly,  it  is  desirable  to  find  an  expression  for  the  phase  function 
and  phase  matrix  in  which  the  dependence  on  scattering  angle  0 is  replaced 
by  dependence  on  incoming  and  outgoing  polar  angles  v and  azimuths  <]> . Thus, 
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if  incoming  and  outgoing  directions  are  denoted  by  (vQ,  <pQ)  and  (v,  <p)  , 
respectively,  they  are  related  to  the  scattering  angle  by 


cos(e)  = yyQ  + (l-y/')1^2(l-y‘-)1/Zcos(<J>-<J>  ), 

KJ  Q STT0 

where  y=cosv  and  y =cosv  • The  corresponding  transformation  for  the 
phase  function  is  obtained  by  expanding  it  in  a Legendre  series: 


1/2 


(A-55) 


N 

p ( e ) = i u£p£(coso), 
£=0 


(A-56) 


with  coefficients  to0  given  by 


»,  - m 

l 2 


p(9)  P (cos0)dcos0 


(A-57) 


and  Legendre  polynomials  P^Ccos©).  When  p(0)  is  normalized,  10  is 
again  the  single  scattering  albedo.  The  transformation  of  p(0^  into 
p(y,  > yQ,  <PQ)  follows  from  the  addition  theorem  for  Legendre  polynomials 
[24]: 


p(y,4sPo>4>0)  - 


N 

> wn 

X, 

£=o 


V“>W 


+ 2 y.,  oSH  p"<i,)p"(i'o)cos",<*-*o) 

m=l 


(A-58) 


Inversion  of  the  order  of  summation  gives 


N 

Z 


N 


p(y,4>,y  , 4>  ) = Z (2-6™)  [ Ew£lt_ElL 
o o o (£+m) ! 

m-o  £=mv  ' 


Pj(W) 


P?(»0  ] cosm ) 

A/  O O 


(A -59) 


where  6™  Is  the  Kronecker  delta.  The  corresponding  expressions  for  the 
phase  mStrix  M are  considerably  more  difficult,  and  their  discussion  will 
be  included  in  a report  which  discusses  multiple  scattering.  The  phase 
function  transformation  (Eqs.  (A-56)  and  (A-57))  is  included  here  since 
it  is  incorporated  in  the  algorithum  used  in  the  ASL  scattering  code. 

The  sequence  of  multipole  moments  , . . . ,b j , . . . is  notoriously  slow 
to  converge  since  at  least  a terms  are  needed.  For  economy  reasons, 
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approximations  are  frequently  used  to  avoid  the  Mie  formalism.  Van  de 
Hulst  [3]  and  Kerker  [4]  give  a number  of  approximations  commonly  used 
and  discuss  their  ranges  of  validity  and  their  error  penalties.  Of  those 
approximations,  Rayleigh  scattering,  valid  for  particles  small  with  respect 
to  incident  wavelength,  is  probably  the  one  most  commonly  used.  Two  forms 
of  this  approximation  are  included  here:  one  which  describes  scattering  by 

small,  discrete  spheres,  and  one  which  describes  scattering  by  gaseous 
molecules . 

Penndorf  [25]  has  derived  approximate  equations  for  the  efficiency  factors 

Q , Q , and  Q by  expanding  the  Mie  multipole  moments  a1 , bj , a2 , 

6xt  s ca  ELD  s 

b , and  a into  polynomials  in  a for  complex  refractive  index  m = n-ik  = 
n(l-ii<) . ^The  result  for  the  extinction  efficiency  factor  is 


^ 24n 2k  . 4n2ic  , 20n2K  , , n ? 

q = a -|-  j — - h — + 4.8n^ic 

ext  z 1 15  3z 


7 (n2+n2ic2)  +4  (n^-n'  -5) 

„ 2 


+ _8  I [ (n2+n2K2)  2+(n2-n2ic*  - 2)  ] 2-36n*  k2  14 
3 \ 2 

1 


1 


>or 


(A-60) 


with 


z = (n'+n'  k2)'  + 4 (n2-n‘ k2 )+4 
1 


(A-61) 


and 

z = 4(n2+n2K2)2  + 12 (n2-n2K2)+9 . 

2 

The  scattering  efficiency  factor  is 

Q = (8/3z2){[  (n2+n2K2)  ‘-hr  -n2K  -2  ] ‘ +36n  V2  la14 
sea  1 


(A-62) 


•{l+(6/5z  ) [ (n2+n2K2)2-4]  or  - (24n2ica3/3z  )}, 
1 1 

and  the  absorption  efficiency  factor  is 


(A-63) 


Q 


abs 
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(A-64) 
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or 


Q . -Im 

abs 


, ,m2-lx 
-4a  ( ) 

m2+2  J 


(A-65) 


For  nonabsorbing  spheres  (k=0),  the  scattering  efficiency  factor  (for 
small  size  parameters)  reduces  to 


= Q, 


sea  ext  3 2.0 

n +z 


8a4  (nizl)2  I ± + 6 2(n!Z2) 


n2+2 


4r  3 ,n6+41n4-284n2+284x 
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(A-66) 


In  quoting  these  approximations,  Kerker  [4]  suggests  that  for  k = 0 
these  equations  may  be  used  for  values  of  a £ 1.4  and  m _<  2.  The  error 
incurred  is  less  that  2%  for  m <_  1.5,  but  may  reach  15%  at  hi  = 2. 

For  absorbing  spheres  the  approximations  are  useful  for  a <■  0 . 8 in  the 
range  1.25  _<  n 1.75  and  nK  <_  1. 

Van  de  Hulst  [3]  has  given  equations  for  cross  sections  for  clouds  of  N 
anisotropic  molecules  with  random  orientations.  For  a gas  with  refractive 
index  m = n-ik,  the  absorption  cross  section  is 


= 47Tk 

ext  AN  ’ 


(A-67) 


which  gives  the  usual  absorption  coefficient  for  gases: 

NC  = 4nk/A.  (A-68) 

ext 

The  corresponding  scattering  cross  section  is  given  by 

C = 8Ti__f — [ (n2-l)  2+4k^  ] , (A-69) 

SCa  3A4N2 
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in  which  f is  the  Cabannes  factor  (f  = 1.054  for  air).  For  k=0,  Eq.  (A-69) 
reduces  to  the  value  given  by  Van  de  Hulst  for  nonabsorbing  gases.  The 
derivation  for  absorbing  gas  cross  sections,  when  done  in  two  apparently 
equivalent  ways,  leads  to  the  4k2  term  being  either  positive  or  negative. 
The  positive  choice  here  is  consistent  with  that  of  Heddle  [26], 


The  theory  given  in  the  previous  portion  of  the  Appendix  has  been  implemented 
in  numerical  form  as  a control  routine,  PGAUSS,  and  several  subroutines,  MIEGSS, 
GAUSS,  VERIFY,  and  GUSSET.  The  PGAUSS  routine  receives  wavelength-dependent 
refractive  index  data  and  definitions  of  particle  size  distributions  from 
the  user,  and  employs  the  subroutines  to  return  scattering  and  extinction 
cross  sections,  phase  functions,  and  phase  matrix  elements  for  the  poly- 
dispersion. PGAUSS  is  an  extremely  versatile  routine  since  it  can  be 
used  to  generate  Mie  functions  for  arbitrary  mono-  and  polydispersions, 
or  a variety  of  standard  types  of  polydispersions  for  comparatively  wide 
ranges  in  particle  size  parameter  and  refractive  index. 

The  PGAUSS  routine  consists  of  three  nested  control  loops.  The  outer 
loop  indexes  GAUSS  quadrature  orders,  the  intermediate  loop  indexes  wave- 
lengths, and  the  inner  loop  integrates  Mie  functions  over  the  size 
distributions  for  each  wavelength.  No  integration  over  wavelength  is 
provided  at  present,  although  it  would  be  simple  to  include  it  as  an 
option. 

The  internal  size  distributions  offered  as  options  include  entries  into 
an  arbitrary,  a log-normal,  the  Deirmendjian  [7]  haze  C,  the  Junge  [6], 
and  the  Khrgian  and  Mazin  [8,  9]  modified  ganma  distributions.  The 
distributions  each  call  for  their  particular  input  data,  and  because  the 
distributions  n(r)  are  normalized  internally, 


the  number  density  information  must  be  input  into  the  model.  The 
resulting  Mie  functions  are  thus  equivalent  to  single  (but  averaged) 
particle  data. 

The  Mie  functions  are  computed  by  the  subroutine  MIEGSS,  which  returns 
extinction  and  scattering  efficiency  factors  Qext  and  Qsca,  intensity 
factors  i^  and  i2,  and  the  relative  phase  6 of  and  i2  for  each 

complex  refractive  index  m and  size  parameter  a.  The  Ricatti-Bessel 
functions  and  their  derivatives  in  Eqs.  (A-5)  and  (A-6)  are  presently 


PGAUSS  PROGRAM 


r 


(A-70) 


r 


1 
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computed  by  forward  recursion  of  the  x?s  and  backward  recursion  of 
the  T?s.  Start  forward  recurring  using  the  following:  [Note:  £ (z)  = 

^n(z)  + ixn(z)] 


XQ(z)  = cos  z. 


(A-71) 


Xj(z)  = cosz/z  + sinz, 

and  the  recursion  formula 

f n ( z ) = (2n-l)  fn_1(z)/z  - fn_2(z) 


and 


fn(z)  “ fn-i(z)  _ (n+1)  fn(z)/z. 


(A-72) 

(A-73) 

(A-74) 


If  L orders  are  needed,  then  the  backward  recursion  is  used  until 
L orders  of  x's  are  generated.  The  recursion  is  continued,  but  now 
the  sum 


S 


z I 

i=n  Xi(z)Xi+l(z> 


(A-75) 


is  formed.  This  is  rapidly  converging  as  soon  as  n is  greater  than  z. 

The  sum  is  converged  when  successive  terms  are  less  than  some  tolerance; 
e.g.,  L = 10  5;  then  from 

(z)  _ 

JL_  Sn>  (A- 

Xn(z) 


we  may  recur  backwards  using  the  relations 

\(z)  = XL(z)‘SL’  (A-77) 

\-l(z)  = xL-l(z)’sT+  > (A-78) 

L Xl^z) 
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and  the  backward  recursion 


Y (z)  = (2n+l)  » (z)/z  - ¥ (z)  (A-79 ) 

n-1  n n4*l 

to  generate  the  V’s.  Finally,  a scaling  factor  F is  calculated  for 

Y (z): 

o 

F = Hf  (z)  /sinz 
o 

and  then  all  the  ¥*s  are  scaled  appropriately.  This  algorithm  is  from 
Wills  [27]  and  appears  to  be  accurate  and  stable  over  the  index  ranges 
of  nka  used  in  this  report. 

The  Mie  series  is  terminated  either  when  two  successive  terms  have  la  1 

1 n 1 

+ |bn|<10“5,  or  when  the  number  of  terms  exceeds  6 4*  Fa.  F is  1.2  for 

a_^51  and  is  1.  4-  2.26a  *613  for  a_>51.  These  cutoff  criteria  terminate 
the  series  satisfactorily  before  serious  instabilities  occur.  The  10“5 
tolerance  should  be  adjusted  to  be  compatible  with  the  computer 
single-word  precision.  This  value  is  reasonable  for  32-bit  machines,  but 
10“k  is  more  appropriate  for  36-bit  machines  and  10-11  should  suffice  for 
60-bit  machines. 

The  angular  functions  tt  (0)  and  t (0)  are  computed  by  forward  recursion 

c n n 

from 

* (9)  = [ (2n-l)  cos0tt  (0)  - nu  (0)]/(n-l),  (A-80) 

n n-i  n-z 


p (0)  = (2n-l)  it  (0)  + p (0), 
n n-i  n-z 


(A-81) 


T (0)  - t\  (0)  cos  0 - sin20p  (0),  (A-82) 

n n n 

and  the  starting  values  tt  = 0,  p = 0,  ir  =1,  and  p =0. 

o o i i 


The  Legendre  expansions  of  the  integrated  phase  functions  and  the  phase 
matrix  elements  are  provided  by  subroutine  GAUSS.  The  integral  in  Eq. 
(A-57)  is  replaced  by  a Gauss -Legendre  sum,  so  that  the  expansion 
coefficients  are  given  by 


21+1 


/ 

-1 

L 


f (cos0)p£ (cos  0)  d cos  0 


2 £4-1  \ 

— 2~  L ^ f (cos0i)P^ (cos0 . )dcos0 ,w. 


i=l 


l li 


(A-83) 


The  replacement  is  exact  (to  machine  accuracy)  when  the  product  in  the 
integrand  is  a polynomial  in  cos  0^  of  degree  2L  - 1 or  less.  The 

order  of  quadrature  L is  an  input  to  PGAUSS  and  should  be  set  at  two  to 
three  times  the  largest  size  parameter  expected  in  the  Mie  calculation. 
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The  abscissa  cos  0^  and  weights  are  computed  by  subroutine  GUSSET. 

The  Mie  angle-dependent  factors  are  computed  at  angles  corresponding 
to  the  cos  0j_. 

In  addition  to  computing  the  coefficients,  GAUSS  reconstructs  p(cos 
as  Pc^c§s  9 and  computes  an  rms  distance  of  p(cos  0^)  from 

Pc(cos  0^)1^  (p (cos  0.)  - p (cos  0,))2]1/2  as  each  successive  term  is 
i=l  1 C 

added  to  the  series  (Eq.  (A-56)).  The  distance  function  is  particularly 
useful  in  determining  the  number  of  terms  to  be  retained  in  Eq.  (A-56)  for 
adequate  precision  in  the  expansion.  Subroutine  VERIFY  produces  a 
reconstruction  pc(cos  0^)  using  the  number  of  terms  which,  minimize  the 
distance  function. 

The  quadrature  weights  and  abscissas  are  computed  by  GUSSET  in  double 
precision  using  the  Davis  and  Rabinowitz  algorithm  [28].  The  abscissas 

= cos  0^n  are  the  k = 1,...,  n zeros  of  Pn(x^n)  = 0,  while  the  weights 
are  given  by 


*kn  " 2d  - *k„>2/["Pn-l<*kn>]2- 


(A-84) 


Initial  estimates  of  the  zeros  are  obtained  from  the  successive  zeros 
of  the  Bessel  function  U0(j^)  = 0 via 


xkn  = cos[jk/((n+l/2)2  + (1  - (2/tt)2)/4)1/2. 


(A-85) 


Final  values  of  the  xkn  are  found  by  Newton-Raphson  iteration. 
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Code  S&E-AERO-Y 
NASA  (Mr.  Vaughan) 

Marshall  Space  Flight  Center, 

Alabama  35802 

NASA  Aero-Astrodynamics  Laboratory 
ATTN:  S&E-AERO-YS 

Marshall  Space  Flight  Center, 

Alabama  35812 

Commander 

USAICS 

ATTN:  ATSLCTD-MS  (Met  Branch) 

Fort  Huachuca,  Arizona  85613 

Commander 

Headquarters,  Fort  Huachuca 
ATTN:  Tech  Ref  Div 

Fort  Huachuca,  Arizona  85613 

Dr.  Robert  Durrenberger 

Director,  The  Laboratory  of  Climatology 

Arizona  State  University 

Tempe,  Arizona  85281 

Commander 

ASL  Yuma  Met  Team 

ATTN:  Captain  Jones 

Yuma  Proving  Ground,  Arizona  85364 

Commander 

Yuma  Proving  Ground 

ATTN:  STEYP-AD,  Tech.  Library 

Yuma,  Arizona  85364 

Fugro  National,  Inc. 

718  East  3rd  Street 
ATTN:  Mr  Jerry  Allen 

Long  Beach,  California  90801 

Sylvania  Elec  Sys.  Western  Div 
ATTN:  Technical  Reports  Library 

P.  O.  Box  205 

Mountain  View,  California  94040 


Geophysics  Officer 
Code  3250 

Pacific  Missile  Range 
Point  Mugu,  California  93042 

Commander 

Naval  Electronics  Laboratory  Center 

ATTN:  Library 

San  Diego,  California  92152 

Director 

Atmospheric  Physics  & Chem  Lab 
Room  31,  NOAA 
Department  of  Commerce 
Boulder,  Colorado  80302 

Library-R-51-Tech  Reports 
Environmental  Research  Labs 
NOAA 

Boulder,  Colorado  80302 

National  Center  for  Atmos  Res 
NCAR  Library 
P.  O.  Box  3000 
Boulder,  Colorado  80302 

E&R  Center 
Bureau  of  Reclamation 
ATTN:  D-751,  Bldg  67 

Denver,  Colorado  80225 

Armament  Dev  & Test  Center 
ADTC  (DLOSL) 

Eglin  AFB,  Florida  32542 

Commandant 

USA  Southeastern  Signal  School 
ATTN:  ATSO-CTD 

Fort  Gordon,  Georgia  30905 

Chief,  Technical  Services  Div 
DCS/Aerospace  Sciences 
ATTN:  AWS/DNTI 

Scott  AFB,  Illinois  62225 


Commander 

US  CDC  Combat  Arms  Dev  Act 

ATTN:  ATCACC-CW 

Fort  Leavenworth,  Kansas  66027 

Commander 

US  Army  Test  & Eval  Cmd 

ATTN:  AMSTE-FA 

Aberdeen  Proving  Ground,  MD  21005 

Commander 

US  Army  Test  & Eval  Cmd 
ATTN:  AMSTE-NBC 

Aberdeen  Proving  Ground,  MD  21005 

Director 

US  Army  Ballistic  Research  Lab 

ATTN:  AMXBR-AM 

Aberdeen  Proving  Ground,  MD  21005 

Commander 

US  Army  Ballistic  Resch  Labs 

ATTN:  AMXBR-IB 

Aberdeen  Proving  Ground,  MD  21005 


Commander 

US  Army  Ballistic  Rsch  Labs 

ATTN:  AMXBR-IC 

Aberdeen  Proving  Ground,  MD  21005 

Director 

USA  Ballistic  Research  Laboratories 

ATTN:  AMXBR-CA 

Aberdeen  Proving  Ground,  MD  21005 

Commander 

US  Army  Armament  Command 
ATTN:  SAREA-DE-N,  C.  G.  Whitacre 

Aberdeen  Proving  Ground,  MD  21010 

Commander 

Aberdeen  Proving  Ground 
ATTN:  STEAP-TL  {Bldg  305) 

Aberdeen  Proving  Ground,  MD  21005 

NASA  Scientific  & Tech  Info  Fac 
ATTN:  Acquisitions  Br 

P O.  Box  33 

College  Park,  Maryland  20740 

Dir  of  Dev  & Engr 
Defense  Systems  Div 
ATTN:  SAREA-DE-DDR 

H.  Tannenbaum 

Edgewood  Arsenal,  APG,  MD  21010 

Commander 
Edgewood  Arsenal 
ATTN:  SAREA-TS-L  (Tech  Lib) 

Edgewood  Arsenal,  MD  21010 

Director  National  Security  Agency 

ATTN.  TDL 

Fort  George  G.  Meade, 

Maryland  20755 

Technical  Processes  Br 
D823 

Room  806,  Libraries  Div  NOAA 
8060  13th  St 

Silver  Spring,  MD  20910 
Commander 

US  Naval  Ordnance  Lab 

ATTN:  Technical  Library 

White  Oak,  Silver  Spring,  MD  20910 

HQ,  ESD(TRI) 

L.  G.  Hanscom  Field 
Bedford,  Massachusetts  01730 

Meteorology  Laboratory 
AFCL(LYV) 

L.  G Hanscom  Field 
Bedford,  Massachusetts  01730 

Air  Force  Cambridge  Rsch  Labs 
ATTN:  LCH  (A.  S.  Carten,  Jr.) 

L.  G Hanscom  Field 
Bedford,  Massachusetts  01730 

Air  Force  Cambridge  Rsch  Labs 
ATTN:  LKI 

L.  G.  Hanscom  Field 
Bedford,  Massachusetts  01730 


US  Army  Liaison  Office 
MIT  Bldg  26,  Rm  131 
77  Massachusetts  Ave. 

Cambridge,  Massachusetts  02139 

US  Army  Liaison  Office 
MIT- Lincoln  Lab,  RM  A-117 
P.  O.  Box  73 

Lexington,  Massachusetts  02173 

The  Environmental  Resch  Institute  of  Michigan 
ATTN:  IRIA  Library 

P.  O.  Box  618 

Ann  Arbor,  Michigan  48107 

Mr.  William  A.  Main 
USDA  Forest  Service 
1407  S.  Harrison  Road 
East  Lansing,  Michigan  48823 

Commander 

US  Army  Tank-Automotive  Cmd 
ATTN:  AMSTA-RHP 

Warren,  Michigan  48090 

Dr.  A.  D.  Belmont 
Research  Division 
P.  O.  Box  1249 
Control  Data  Corp. 

Minneapolis,  Minn.  55440 

Director 

USA  Engr  Waterways  Exper  Sta 
ATTN:  Library  Branch 

Vicksburg,  Mississippi  39180 

Director 

LISAE  Waterways  Experiment  Station 
ATTN:  Mr  A.  N.  Williamson 

P.  O.  Box  631 

Vicksburg,  Mississippi  39180 
Commander 

US  Army  Cold  Regions  R&E  Lab 
ATTN:  CRREL-RP 

Hanover,  New  Hampshire  03755 

Commanding  Officer 
Picatinny  Arsenal 
ATTN:  SARPA-TS-S,  #59 

Dover,  NJ  07801 


Commander 

US  Army  Satellite  Comm  Age 
ATTN  AMCPM-SC-3 
Fort  Monmouth,  NJ  07703 


Commander 

US  Army  Electronics  Command 
ATTN-  AMSELBL-CT-S 
Fort  Monmouth,  NJ  07703 


Commander 

US  Army  Electronics  Command 
ATTN  AMSEL-CT-D 
Fort  Monmouth,  NJ  07703 


Commander 

US  Army  Electronics  Command 
ATTN  AMSEL-GG  TD 
Fort  Monmouth,  NJ  07703 


Commander 

US  Army  Electronics  Command 
ATTN*  AMSEL-MS-TI 
Fort  Monmouth,  NJ  07703 

Commande  r/ Di  rector 
US  Army  Electronics  Command 
ATTN:  AMSEL-NL-D 

Fort  Monmouth,  NJ  07703 

Commander 

US  Army  Electronics  Command 

ATTN:  AMSEL-NL-H  (EVANS  AREA) 

Fort  Monmouth,  NJ  07703 

Commander 

US  Army  Electronics  Command 
ATTN:  AMSEL-RD-D 

Fort  Monmouth,  NJ  07703 

Mr.  Edward  C.  Nell 
US  Army  Electronics  Command 
ATTN:  AMSEL-SI-CB 

Fort  Monmouth,  NJ  07703 

Director,  Electronics  Technology  and  Devices  Lab 
US  Army  Electronics  Command 
ATTN:  AMSEL-TL-D,  Bldg  2700 

Fort  Monmouth,  NJ  07703 

Commander 

US  Army  Electronics  Command 
ATTN:  AMSEL-VL-D 

Fort  Monmouth,  NJ  07703 

Commander 

US  Army  Electronics  Command 
ATTN:  AMSEL-WL-DI 

Fo/t  Monmouth,  NJ  07703 

^ir  Force  Weapons  Lab 
ATTN:  Technical  Library  (SUL) 

Kirtland  AFB,  NM  87117 

Commander 

White  Sands  Missile  Range 
ATTN:  Technical  Library 

White  Sands  Missile  Range, 

New  Mexico  88002 

CO,  US  Army  Tropic  Test  Center 
ATTN:  STETC-MO-A  (Tech  Lib) 

APO  New  York  09827 

Rome  Air  Development  Center 
ATTN : Documents  Library 

TILD  (Bette  Smith) 

Griffiss  Air  Force  Base, 

New  York  13440 

US  Army  Rsch,  Off.  Durham 
Box  CM,  Duke  Station 
ATTN:  CRDARD-IP 

Durham,  North  Carolina  27706 

Environmental  Protection  Agency 
Meteorology  Laboratory 
Research  Triangle  Park, 

North  Carolina  27711 

Commandant 

US  Army  Field  Artillery  School 
ATTN:  ATSF-CTD-A 

Fort  Sill,  Oklahoma  73503 


Commandant 

US  Army  Field  Artillery  School 
ATTN:  Met  Division 

Fort  Sill,  Oklahoma  73503 

President 

US  Army  Field  Artillery  Board 
Fort  Sill,  Oklahoma  73503 

Field  Artillery  Consultants 
1112  Becontree  Drive 
ATTN:  COL  Buntyn 

Lawton,  Oklahoma  73501 

Air  Force  Avionics  Lab 
ATTN:  AFAL/TSR 

Wright-Patterson  AFB,  Ohio  45433 

Commander 
Frankford  Arsenal 
ATTN:  Library,  K2400,  Bldg  51-2 

Philadelphia,  PA  19137 

Commander 
Frankford  Arsenal 

ATTN  SARFA-FCD-O,  Bldg  201-2 
Bridge  & Tarcony  Sts. 

Philadelphia,  PA  19137 

Commander 
Frankford  Arsenal 
ATTN:  W 1000-65-1 

Philadelphia,  PA  19137 

Commander 

US  Army  Nuclear  Agency 
Building  12 

Fort  Bliss,  Texas  79916 
Commander 

US  Army  Air  Defense  School 
ATTN:  C&S  Dept,  MSLSCI  Div 

Fort  Bliss,  Texas  79916 

Commander 

US  Army  Dugway  Proving  Ground 
ATTN:  STEPD-AD-L  (Tech  Library) 

Building  5330 
Dugway,  Utah  84022 

Commander 

US  Army  Dugway  Proving  Ground 
ATTN:  TO-M 

Dugway,  Utah  84022 

Mr.  Frank  D.  Eaton 

Department  of  Soil  Science  and  Biometeorology 
Utah  State  University 
Logan,  Utah  84322 

Commander 

US  Army  Materiel  Command 
ATTN:  AMCRD-OS  (Mr.  Andrew) 

Alexandria,  Virginia  22304 

Mr.  Frederick  L.  Homing 
5001  Eisenhower  Avenue 
ATTN:  AMCRD-TV 

Alexandria,  Virginia  22304 

Commander 

USA  Security  Agency  Combat  Dev  Actv 
Arlington  Hall  Station,  Bldg  420 
Arlington,  Virginia  22030 


Commander 

US  Army  Security  Agency 
ATTN:  IARD-T 

Arlington  Hall  Station 
Arlington,  Virginia  22030 

CO,  USA  Foreign  Sci  & Tech  Center 
ATTN:  AMXST-ISI 

220  7th  Street,  NE 
Charlottesville,  Virginia  22901 

Commander 

US  Army  Engineer  Topographic  Lab 
(STINFO  CENTER) 

Fort  Belvoir,  Virginia  22060 

Commander 
Eustis  Directorate 
US  Army  Air  Mobility  R&D  Lab 
ATTN:  Technical  Library 

Fort  Eustis,  Virginia  23604 

Commander 

US  Army  Combat  Dev  Cmd 
ATTN:  ATCD-CI-I 

Fort  Monroe,  Virginia  23651 

Commander 

HQ  US  Army  Training  & Doctrine  Cmd 

ATTN:  ATCD-FA 

Fort  Monroe,  Virginia  23651 

Marine  Corps  Dev  & Educ  Cmd 
Development  Center 

ATTN:  Cmd,  Control,  & Comm  Div  (C>) 

Hochmuth  Hall 
Quantico,  Virginia  22134 

Commander 

US  Army  Arctic  Test  Center 
ATTN:  STEAC-PL-TS 

APO  Seattle  98733 

Director,  Systems  R&D  Service 
Federal  Aviation  Administration 
ATTN.  ARD-54 
2100  Second  Street,  S.W. 

Washington,  DC  20590 

Chief,  Weather  Climate 
NASA  Headquarters 
ATTN  Code  ERM 
Washington,  DC  20546 

Director 

Naval  Research  Laboratory 
CODE  2627 
Washington,  DC  20375 

Dir,  US  Naval  Research  Lab 
ATTN : Code  5530 

Washington,  DC  20375 

Commander 

Naval  Elect  Sys  Cmd  HQ 
ATTN:  CODE  51014 

Washington,  DC  20360 

The  Library  of  Congress 
ATTN:  Exchange  & Gift  Div 

Washington,  DC  20540 


Commanding  Officer 
Harry  Diamond  Laboratory 
ATTN:  Library 

Washington,  DC  20438 

National  Weather  Service 
National  Meteorological  Center 
World  Weather  Bldg  - 5200  Auth  Rd 
ATTN:  Mr  Quiroz 

Washington,  DC  20233 

Director 

Defense  Nuclear  Agency 
ATTN:  Tech  Library 

Washington,  DC  20305 

David  H Slade 

Leader,  Physical  and  Chemical  Transport  Pgrm 
Division  of  Biomedical  and  Environmental  Rsch 
ATTN:  US  Atomic  Energy  Commission 

Washington,  DC  20545 

Commander 

Naval  Air  Systems  Command 
Meteorological  Division 
AIR-540 

Washington,  DC  20360 
Commander 

Naval  Weather  Service  Command 
Washington  Navy  Yard 
Bldg  200,  Code  304 
Washington,  DC  20374 

Commander,  Naval  Ship  Systems  Cmd 
Technical  Library,  Rm  3 S-08 
National  Center  No.  3 
Washington,  DC  20360 

Dr  John  L.  Walsh 
Code  5503 
Navy  Research  Lab 
Washington,  DC  20375 

Head,  Atmospheric  Sci  Section 
National  Science  Foundation 
1800  G.  Street,  N.W. 

Washington,  DC  20550 

Head,  Met  Rsch  Branch 
Meteorological  Division 
Naval  Air  Sys  Cmd  (AIR-5401) 

Washington,  DC  20361 

HQDA  (DAEN  RDM/Dr.  de  Percin) 
FORRESTAL  BLDG 
Washington,  DC  20314 

Ofc  Asst.  Sec  of  the  Army 
R&D  ATTN : Asst,  for  Rsch 

Rm  3-E-379,  The  Pentagon 
Washington,  DC  20310 

Military  Assistant,  Environmental  Sciences 

ODDR&E/E&LS 

Pentagon,  Room  3D  129 

Washington,  DC  20301 

USAFETAC/CB 

Navy  Yard  Annex,  Bldg  159 

Washington,  DC  20333 


Air  Force  Weapons  Lab 
ATTN:  Technical  Library  (SUL) 

Kirtland  AFB,  NM  87117 
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